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ABSTRACT 


In one method of solution to the problem of making up a schedule 
of work in a job-shop" type industry there arises a symmetric matrix--- 
the conflicts among the jobs to be done, This thesis is primarily 
concerned with the transformation of this "conflict" matrix to another 
matrix---a diagonalized symmetric matrix, That is a matrix having its 
nonzero elements bunched or clustered near the principal diagonal, 
The transformation is considered as the product of simple transformations 
each of which interchange rows or columns of the undiagonalized or 
partially diagonalized matrix, 

A numerical technique is devised for the systemmatic reduction 
of an undiagonalized symmetric matrix to successively better diagonalized 
matrices, The technique is programmed for use on the NCR 102A Digital. 


Computer, 


T 
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CHAPTER I 
INTRODUCTION TO THE PROBLEM UNDER CONSIDERATION 
l. Introduction 

In a "job-shop" type industry there are many jobs---job orders--- 
to be Jone, each requiring the use of certain of the varied facilities, 
special tools, and specially trained technicians available in the 
industry, The question of how to make up the time schedule of work 
to be done arises, That is which jobs should be time scheduled first 
considering the complexity and use of varied facilities, etc,, in the 
industry, Then having time scheduled the first one, which job should 
be time scheduled next, and so forth, 

The mathematical problem to be considered in this thesis will assist 
in the determination of an orderly procedure for time scheduling the jobs 
to be done and will provide the order in which jobs should be considered 
for time scheduling, 

2. The Mathematical Problem 

Given a symmetric matrix A, consisting of zero and nonzero elements; 
find a symmetric matrix С whose elements satisfy the following; 

1) The elements of a given column of A become the elements of 
some column of C, 

11) The elements of a given row of A become the elements of some 
rows of C, 

11i) The nonzero elements of the matrix C are located along the 
principal diagonal of the matrix, That is to say there is a clustering 


of elements near the principal diagonal, 





iv) In addition be able to enumerate the correspondence of columns 
of the matrix A with columns of the matrix C, 

The conditions i and ii together are equivalent to the 
following: 

i!) Columns may be rearranged, maintaining the order of elements, 

11!) Rows may be rearranged, maintaining the order of elements, 

For consider any element of the matrix say (i, j). The column i is 
shifted to say the h-th column not necessarily retaining the order of 
elements, Also row j is shifted to say the k-th row not necessarily 
retaining the order of elements, The element that was in position 

(1, 3) must now be in position (h, k), Houever, the same result is 
obtained by shifting the i-th column to the h-th column maintaining the 
order of elements, and then shifting the j-th row to the k-th row main- 
taining the order of elements, This must be true for all (1,j). Thus 
i and ii imply i' and ili', That і! апа їі! imply i and ii is obvious, 

Àn assumption is made that the number of nonzero elements in any 
row or column is small compared to the order of the matrix, say one tenth 
the elements are nonzero, If this is not the case the idea of "diagonal- 
ization" or clustering the elements along the principal diagonal has little 
significance, 

The transformation or sequence of transformation used to find the 
matrix C should be applicable to matrices of large order say 62 x €2 or 
even as large as 500 x 500, Since we are considering large order matrices 
it seems logical that the transformation or sequences of transformation 
should be susceptible to programming on a computer, say the NCR 102A 


Digital Computer, 





3, Notation, 
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Figure 1 General Form of Matrix 


A symmetric matrix A of order n, ann x n array of zero and non- 
zero elements, The exact value of the nonzero elements 1s not important-- 
the reasons for this are explained in Chapter IV, 

All nonzero elements Ру; are located by coordinates (Cy »R,) 
indicating its position in the matrix array is in the C, th column 


and the R,th row, All elements not designated and enumerated are 


j 


assumed zero, 


The columns are designated 07,6 enn. numbered in the usual 


As 
fashion from left to right, 





The rows are designated ОО: numbered from the bottom 
upward, This unusual numbering is associated with the numerical technique 
used to "diagonalize" the matrix. 

The sum of the row and column positions is designated 514* where 
533 EG Rs. Note that all the elements of the principal diagonal 
have the same value, namely, n * 1, 

The distance, D, of an element from the principal diagonal is 
E -(n + 1)| =D, This is an integer distance along the ith column or 
the jth row, 

4. The Measure of Diagonalization öf a Matrix 

Some criteria must be set up to measure the clustering or bunching 
of the nonzero elements along the principal diagonal to distinguish 
between the matrices, That is to determine which of two matrices is 
more diagonalized, This measure is relative and serves only as an 
index, 

A simple criterion would be to minimize the distance of the element 


farthest from the principal diagonal, 


С (п + 1) “Ss, | -——-Min Max D 


This will produce a certain amount of diagonalization , However, a 
simple example will show that a more powerful criterion is needed, 
Figures 2a and ?b both have the same Min Max D, namely 2, however, 


Figure 2b is more compact along the principal diagonal, 
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Figure 2 Examples of Min Max D Criteria 
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Consider the criterion that the sum of the distance from each 


element to the principal diagonal squared be minimized., 


Min » [tn De з, ў -—-Mn Ў р” 
Р,, 
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Each element far from the principal diagonal will contribute greatly 
to the summation, whereis elements next to the principal diagonal will 
have negligible effect on the summation, Thus the greater the cluster- 
ing along the principal diagonal, the smaller will be the summation, 
The previously considered Min Max D criterion will tend to be satisfied 
also, For consider Figure 3 showing two columns (all other columns are 
assumed the same in both cases and, therefore, contribute the same to 


the summation in both cases), 
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Figure 3 Example of Г р” for Row Interchange 


In the summation note that only the contribution of the outer element 
in each column is affected by the interchange of the columns under con- 
sideration in Figure 3. In (a) the outer elements contribute to the 
summation (k«1)^ * (1-1) or 2242. In (b) these same elements con- 
tribute 2k^. Thus (b) is better diagonalized in the sense of both 

Min Max D and Min 7 D^, 





For the remainder of the discussion and in testing the merits of 


a transformation under consideration we will use Min = 








CHAPTER II 


SOME FEATURES OF MATRICES AND ELEMENTARY TRANSFORMATIONS 
A square matrix A whose elements satisfy the conditions 


a = э 


1] 11 
possesses symmetry about its principal diagonal, and is called a 
symmetric matrix. 

15% C=BA, then the matrix A is said to be linearly trans- 
formed into the matrix C, The matrix B, which affects the transformation 
of A, may sometimes be advantageously thought of as the multiple product 


of very simple matrices, 
B BEE ¿BED 


The transformation of A is thus accomplished through a succession of 
simpler transformations, the first of which is В), the second Bas and 
ӨЫ forth, 

The matrix can be transformed through postmultiplication as well as 
by premultiplication, 1,e., C=AB', In which case the matrix 


affecting the transformation is decomposed into components 


В' = В В,.........В. 


1 


1 
Where B. indicates the first simple transformation, B, 


the second, and 
so forth, 

It is of interest to determine the simplest fundamental forms into 
which can arbitrary nonsingular matrix B may be decomposed, and to 


interpret the individual transformations which they effect upon the form 





of the matrix A, These are called elementary transformations and the 
matrices which produce them are called elementary transformation matrices, 
There are three fundamental types of elementary transformations and 
correspondingly three fundamental types of transformation matrices, 

The first type amounts to an interchange of two rows or columns of the 
matrix A. The second type is the addition of the elements of a row or 
column to the corresronding elements of another row or column, The third 
type is the multiplication of any row or column of A by an arbitrary 
nonzero factor. For our purpose it is necessary to develop only type 

one to a greater extent, 

The transformation may be affected by multiplying A by a transforma- 
tion matrix. If the desired transformation is intended to affect the 
ROWS, A is PREmultiplied by the transformation matrix; if the COLUMNS 
are to be affected, A is POSTmultiplied by the transformation matrix, 

The elementary transformation matrix is formed from the unit matrix 
by performing on it the same elementary transformation which the desired 
elementary transformation matrix is intended to effect in the matrix A by 
means of (pre-)postmultiplication, Thus I (the unit matrix) with its 
pth and qth rows (columns) interchinged yields a transformation matrix 
which by means of pre- (post-‘multiplication interchanges the pth and qth 
rows (columns) of the matrix A, Apparently the transformation matrix 
has two forms depending on whether it is intended to effect a transforma- 
tion of the rows or columns of A, however, it is to be noted that one 
form is the transpose of the other, and further for type one the two forms 
are identical, «е will denote this elementary transformation matrix as 


E where the subscripts indicate the rows (columns) being interchanged, 


8 








As an example, for a fourth order matrix to interchange the second and 


third rows (columns), 


| 

| "23 23 32 

| Columns 

Consider a symmetric matrix A, pre- and postmultiplied by the 


elementary matrix T The resulting matrix 


i-j' 


is a symmetric matrix with both the ith and jth rows and columns inter- 
changed, Continuing this process we may by the elementary matrix 


Tip Produce the matrix 


Ten en Дош ا و‎ 


Thus by a multiple product of elementary transformation matrices pre- 
and postmultiplying the matrix A, the matrix A is linearly transformed 
into a symmetric matrix C, 

The matrix C which may be obtained from the matrix A bv means of a 
finite number of elementary transformations is said to be equivalent to A, 
The equivalence of matrices is a mutual relationship since the trans- 
formations are nonsingular and hence reversible, Thus if C can be 
obtained from A by a succession of elementary transformations, it 
follows that A can be obtained by means of elementary transformations 
from C, The state of equivalence of the two matrices A and C may be 


stated 





PAQ=C where P and Q are nonsingular matrices, 
"ow if A is a symmetric matrix ani Q - Q2۰...” where each 


% is an elementary transformation matrix T, if further 


Ра З а...) БИШ 0, 


then С is also a symmetric matrix, Since 
С = gi AQ 


The transpose of C is then 


T eT 


ch = qi at (qi! 


= Q AQ. 


A matrix which is its own transpose is symmetric, So C is symmetric, 


The preceding discussion has covered the elementary transformation 
matrices in general. Now to apply these concepts in the notation to be 
used in the numerical technique to be developed in Chapter III, 

In order to have a symmetric matrix after premultiplication by 
an elementary transformation matrix P which interchanged the ith and 
jth columns we must multiply by an appropriate row interchanging elemen- 
tary matrix transformation, Since the rows are designated from the 
bottom up, we shall designate the appropriate postmultiplying elementary 


matrix as po Where 
i 


{+1 (п + 1) 


j + 3 = (а +1) 


10 








Examination of Figure 1 will show that this indeed is the appropriate 
row interchaneing matrix, 

Thus A, = C A) 777. 
As an example consider the matrix A interchanging first the 2nd and 


5th columns, then the corresponding rows, 


10101 1 1100 RICO O 
01101 01101 201.001 
A=11110| A=10111 E 10111 
100110 100110; 00110 
OA 11001) 01101 
psc 
AU cM 


ten i =6-2=4 
j=6-5=1 
у ‚лыт, 
A = (T À) Ti 


Similarly it can be seen that if the postmultiplication elementary 
transformation matrix is performed first then followed by the appropriate 


premultiplication elementary transformation matrix we have 


A =T (АТ 


27 tij ) 


ij 


As an example consider the matrix A interchanging first the 2nd and 


5th rows, then the corresponding columns, 





1 0101 Bro 9 1010 o 
01101 Ono DEM 0 l 
А=11110) at =|12210 A, =|1 111 o! 
00110 10101 0011 1 
а 100 1 moo 01011 


PT 





25 
= jo m ! 
A, T A Ty A 
= ) 
Ay = 1,5. | 4A Toe. 
Note that T applied as 4 premultiplication matrix and as a post- 


multiplication matrix produces different transformed matrices, That is, 


Т1 (А75 ) = ag FA = (Т,5А) 7; 


Since we are discussing transformations of symmetric matrices only, 
necessirily these elementary matrix transformations will occur in groups 
of two, So for convenience define a pair of elementary matrix transforma- 
tions as the transformations Ti; applied to the rows(columns) and the 


corresponding Tij applied to the columns (rows), 
The given matrix A has a certain 5 р“, say (2. UR which in the 


previous example was 62, If we apply the transformations 


A, = (Ts, A) Тү; 


the transformed matrix will have a new > De which for the previous 


example is 32, If on the other hand we apply the transformations 


Ay Tr (A Tij) 


a different transformed matrix is obtained with a certain = р? which 
for the previous example is 32, Thus it can be seen that different pairs 
of elementary transformation matrices will produce transformed matrices 
"diagonal ized" to a greater or lesser degree, Some pairs of elementary 
matrix transformations may even produce a transformed matrix which 1s 


less diagonalized in the sense of Min Da 


eun 


12 


CHAPTER III 


PROCEDURE 


1, General 

Basically the procedure is to apply a pair of elementary trans- 
formation matrices to the undiagonalized or partially diagonalized matrix 
and test the diagonalization in the sense of minimizing 2m. 

a, If the transformed matrix is better diagonalized, then use this 
matrix as the matrix to be diagonalized and repeat the process applying 
elementary transformation matrices to this partially diagonalized matrix, 

b, If the transformed matrix is NOT better diagonalized, apply a 
different pair of elementary transformation matrices to the undiagonalized 
matrix and test for Er again, 

We can choose the pair of elementary transformation matrices either 
by some random process or systematically, This thesis is concerned with 
a systematic procedure for improving the diagonalization of a given 
matrix, 

2. A systematic Procedure 


Án element, let us designate it P which lies farthest from the 


IJ’ 
principal diagonal individually contributes the most to > р? of any 
element Pay in the matrix, Since by symmetry there are two such elements 
let us consider the one with the largest S, that is 517° It appears that 
a pair of elementary matrix transformations which involve this element 
will be a good transformation to hypothesize making, Furthermore, to move 
this element to a position one unit from the principal diagonal will 
reduce the contribution tə =: of this element to a minimum, However, 


other elements of the rows and columns interchanged will change distance 


13 





from the principal diagonal so that one must test this hypothesized trans- 
formation to see if it improves the diagonalization, If it does we make 
the hypothesized transformation and repeat the process again choosing 
the new Pry with a maximum 517. 
If the previously hypothesized transformation does not improve the 
diagonalization we hypothesize another transformation, If the previous 
elementary transformation pair TiJ was a rremultinplication first, i.e., 
row interchange first, we now hypothesize i making & column interchange 
first such that Pry is moved to one unit from the principal diagonal, If 
this results in a better diagonalized matrix we make this hypothesized 
transformation and repeat the process again, choosing the new Pry with a 
maximum 517; 
If the hypothesized transformation moving Pry to a position one 
unit from the principal diagonal does not improve the diagonalization of 
the matrix, hypothesize moving Pry to two units from the principal 
diagonal, If neither of these produce a better diagonalized matrix, 
hypothesize moving Prs to three units from the principal diagonal and 
continue until the interchange is with the row (column) adjacent to the 
element Poy: 


If all these hypothesized transformations with this element P., fail 


IJ 
to produce & better diagonalized matrix disregard this element and choose 
from among the remaining elements the element with the largest 5,4. Repeat 
the process of hypothesizing transformations with this element until one 
is found which improves the diagonalization, If none is found among all the 


elements with EM greater than n + 1, i,e., we have looked at all possible 


moves of all ele ents on one half of the symmetric matrix, then the matrix 


14 





which we are attempting to diagonalize is itself the best possible 
diagonalization in the sense of minimizing P 

By hypothesizing moving this "wildest" element to one unit from 
the principal diagonal first, we are maximizing the possible improvement 
in ons stage of diagonalization, The actual improvement achieved іп опе 
stage will usually be somewhat less however, The distance between the 
rows (columns) being interchanged must be variable since it can be shown 
that a system which hypothesizés making the interchange always with an 
adjacent row (column) will not achieve maximum diagonalization, 
Figure 4a is the best diagonalization possible using a hypothesis rule 
of always making the interchanges with the adjacent column (row), 
Figure 4b is the best possible diagonalization using the systematic 


procedure outlined above on the same matrix, 


1⁄5 10000000 E 100000000 
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Figure 4 Example of Systemmatic Procedure For 
Better Diagonalization 


3. Systematic Procedure Applied to the NCR 102A Computer 
All column and row numbering is done in the octal number system 
for the simple reason that the NCR 102A computes in binary number 


system using data inserted in the octal number system, This 


T5 
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information could be inserted as decimal numbers, However, before 

the computer can start to compute a subroutine would have to be 

employed to convert all data to octal, 

The form and amount of information about a particular nonzero 
element of the matrix A which will be placed in the computer is dictated 
in part by the computer and in part a desire by the writer to contain in one 
computer storage word all information relative to this particular element, 
Since a computer storage word is composed of 12 octal digits it is conven- 
ient to store the information about a particular Pj as follows; 


i) for matrices of order 63 x 63 or less 





e. Rio ог; С. К, 


Since numbers as large as 63 сап be represented in 2 octal numbers (63 
decimal is 77 octal) this representation is possible, Where 


in the matrix before 


C Column designation of the nonzero position Pai 


Lo 
any diagonalization has been effected, 


R Row designation of the nonzero position P,, in the matrix before 


jo ij 


any diagonalization has been effected, 
4 Current column designation of this Pag is the matrix at any stage 
of diagonalization, 


j Current row designation of this Р, in the matrix at any stage 


J 
of diagonalization. 


ЕТ) The current С, + К. 


It is to be noted that the C ve and 5, may change from time to time 


1” j 


during the diagonalization process, However, the C 


T and Rio remain 


17 


unchanged and thus a one to one correspondence between the rows (colums) 
of the original matrix and the current partially or completely diagonal- 


ized matrix is obtainable by comparing the R, and R, (C 


Jo J lo 
The detailed numerical techniques described later are based on 


and C,). 


stored word structures for order 63 x 63 or less, 


ii) for matrices of order 511 x 511 or less 





Numbers as large as 511 can be represented in 3 octal numbers (511 decimal) 
is 777 octal), However, certain information must be sacrificed in order 
to contain the essential information in one computer storage word length, 
For these larger matrices the value of Э; must be formed wherever needed 
rather than being immediately available, The symbols have the same 


meaning as above, 
Determination of the largest Si 
In order to determine the "wildest" element or more precisely the 


element with the largest Sy it is necessary to order the elements Р,, 
or at least order some of them, say the largest 6, The reason that more 
than one largest 914 must be obtained is that it may happen that none of 
the hypothesized transformations using the Pry element with the largest 


S++ will improve the diagonalization, However, it is unlikely that 


IJ 


among the elements with the 6 largest 394,2 at least one will not have a 


transformation which will improve the diagonalization, If one of them 


does improve the diagonalization the resulting P,,'s will be reordered 


ij 
18 
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and the 6 largest S,.'s of the new matrix obtained, It is to be noted 


і) 
that the techniaue of ordering tends to order not only the 6 largest but 
to relatively orier all the elements. Thus after several partially 

diagonalized matrices are obtained, all the ELS will be relatively well 


ordered, 


The technique used is simply to compare the 511 of the K and K + 1 


J 
of the K + 1 and K + 2 storage cells, putting in 


storage cells, putting the Psy with the smaller 5; in the K storage cell, 


Then compare the Sj 
K + 1 the n with the smaller je Continuing this process until, the 


N (last) storage cell is reached, then put both Р, у'5 away. Decrease 
the number of storage cells under consideration by one, that is, for the 
second time through this ordering process only the first N-1 storage cells 


wouli have their SE compared. To determine if the preassigned number 


1j 5 have been obtained, test if N-6 is the last storage cell 


to be considered for the next time through the ordering process, If not 


of largest S 


repeat the process starting with the first two storage cells; if so go 


on to the next segment of the program, 


Making a Hypothesized Transformation 
For the first hypothesized transformation we consider interchanging 


the elements of two rows first, The Pry stored in the N storage cell will 


have the largest S It is easily determined which row this P}, ís in, 


IE IJ 
namely the Jth row, and the distance from the principal diagonal, namely D. 
Then J-D is the row which has the Ith column on the principal diagonal, 

Let K=J-D+l 


Then the interchange of the Jth and Kth rows will place P,, next to the 


IJ 
principal diagonal, Now to maintain symmetry two columns must also be 


interchanged, It is easily seen that the columns 
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Figure 7 
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сз! 
il 


(n*1)-J 


(n*1)-K 


хі 
11 


are the appropriate columns to interchange to maintain the symmetry. 
These interchanges are the equivalent of the elementary matrix trans- 


formation 


TI AT IK. 


Forming these quantities J, K, J, K where they are readily available 


proceed to examine each а performing the following tests; 


1) if ] 2: 7 then make j=K 


if j = K then make j = J 
otherwise row designator remains unchanged. 


ii) if i =J then make і = К 


il 
Cal 


if i =K then make i 


otherwise column designator remains unchanged, 


2 


Determine the new S,, of this element and then D, Form D” and add the 


1j 
contribution of this element to the partial zx p* already obtained, 


When all the P,,'s have been so treated, test to determine if the 


1] 
IH p^ for the hypothesized change is less than 5 р“ for the matrix 
before transformation, If indeed this new = De is less, we have a 
matrix more diagonalized than before and we start the process over again, 


ordering by the S,,'s, etc, 


ij 
If, however, the new 2 р is not less than before discard this 
hypothesized transformation and replace the storage cells with the 
partially рае set of PINE which had previously been stored separately. 
This time hypothesize a transformation by interchanging elements of two 


columns first, Taking again Pry in the N storage cell determine which 
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column Pry is in, namely, I, and the distance from the principal diagonal, 
Then I - D is the column which has the Jth row on the principal diagonal, 
Let L-I = D+] 


Then the interchange of the Ith and Lth columns will place P., next to 


IJ 
the principal diagonal, Now to maintain symmetry two rows must also be 


interchanged, It is easily seen that the rows 


I=(n+1)-I 


t~i 
П 


(n +1) -L 


are the appropriate rows to interchange to maintain the symmetry. These 
interchanges are equivalent to the elementary matrix transformation 

pair Tit А тн. | 
Forming these quantities I, 1, T L where they are readily available 


proceed to examine each P,, performing the following tests; 


1j 
i) if j = I then make j =L 
if j = L then make j = 1 


otherwise row designator remains unchanged, 
11) ifi=I then make i=L 


ү 


if i =L tren make i 


otherwise column designator remains unchanged., 


2 


Determine the new 5, of this element and then D, Form D” and add the 


J 
contribution of tris element to the partial T p already obtained, 


When 311 the Р, have been so treated test to determine if the 


j 
E^ р” for the hypothesized change is less than Z р” for the matrix 
before transformation, If indeed this new RT p^ is less, we have a matrix 


more diagonalized than before and we start the process over again, ordering 
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bv the S,4's, etc, 


1j 

If, however, the new È p? is not less than the = р” for the matrix 
before transformation, discard the hypothesized transformation, replace 
the storage cells with the partially ordered set of Pay? which had 
previously been stored away, 
Now make K' = K + 1 

L' =L+]1 

ІР К! < J hypothesize making interchanges as before using K' and L' 
in place of K and L, 

If, however, K' = J, it means that we have hypothesized all possible 
row (column) interchanges which shift the Jth row (Ith column) and none 
of them improve the diagonalization of the matrix, Make N'-N-1, 
that is repeat the hypothesizing transformations process using the next 
wildest element, If none of these hypothesized transformations will 
improve the diagonalization of the matrix, make N' = N - 2, and so on, 
continuing this process until all elements have been tried. If none of 
the hypothesized transformations improve the diagonalization of the matrix, 
then the best diagonalization possible using the criteria of minimize 


EE has been attained. 








CHAPTER IV 


A PHYSICAL PROBLEM 


Consider & job-shop type industry such as a shipyard which may be 
considered typical, The products may be considered the accomplishment 
of job orders, each of which may require the use of special facilities, 
special machines and equipment, and the special talents of specialized 
technicians, 

Let us consider several typical or probable job orders: 

Job Order #1: Main Battery Range finder; clean, repaint pris- 
matical marks on lenses, reevacuate rangefinder tube, fill evacuated 
tube with nitrogen, remount and realign the rangefinder, 

Job Order /2: Main Battery Fire Control Equipment; correct 
systemic errors in synchro system between radar and optical portions 
of fire control system, also accomplish electronic modification XX. 

Job Order 43: Binoculars; collimate and adjust 60 pairs of 
binoculars, 

Job Order #4: Radio Equipment, overhaul model YY receiver, 
fnd so forth...» 

The job orders are analyzed with respect to special machines, 
facilities, and ee technical work required. For example, Job 
Order #1 will require, 1) all the facilities of the optical section for 
work space during cleaning, 2) a specialized technician to repaint the 
prismatical marks on the lenses, 3) special equipment to evacuate the range 
finder tube afterwards, 4) special equipment to refill the range finder 
tube with nitrogen, 5) special alignment equipment and technical assistance 


during remounting, 
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Job order #2, while primarily the concern of electronics technicians 
and requiring electronics equipment also will require a synchro expert 
and technicians from the optical shop, 

Job order #3 requires the shop facilities of the optical shop, Thus 
each job order can be analyzed and a chart prepared as indicated in 
Figure &, where there are m job orders and n columnar headings, 

Figure 8 in effect represents a m x n matrix of zero and nonzero 
elements, Each x in the figure represents a nonzero element, All the 
nonzero elements of a given rov represent the columns (men, machines and 
facilities) which are affected by a particular job order, Post multipli- 
cation of thic m x n matrix M with an n x m matrix, say МІ - the matrix 
transpose- will produce an m xm matrix A, The matrix multiplication ís 
carried out in the usual way using the conventions 


(nonzero element) x (nonzero element) = (nonzero element) 


(nonzero element) x ( zero element) = ( zero element), 

This matrix M M is symmetric with the rows and columns representing 
job order numbers, Thus a nonzero element in the matrix position (i,j) 
indicates that job order i ind job order j have some conflict between them, 
This conflict may be that both job orders require the use of the space 
facilities of the optical shop or the use of the special nitrogen filling 
equipment, This matrix product M M! is thus the "Job order Conflicts matrix, 

This symmetrica! conflict matrix is the matrix which is diagonalized 
elsewhere in this thesis, The diagonalization provides an ordering of the 
job orders such that all job orders which conflict with a given job order 
are listed in the ordering near that given job order, Thus a job order 


which conflicts with no other job orders will appear at one end of the ordering. 
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Note that the ordering is not unique for the ordering may be enumerated 
in reverse order for the same diagonalized matrix, 

Let us consider again the original job order vs men, machines and 
special facilities matrix M, Since there is nothing special about 
the order that the job orders are listed, let us rearrange the rows 
so that the job orders which conflict with each other are near each 
other, That is precisely the ordering of the job orders given by the 
diagonalization of the "Job Order conflict" matrix, When the matrix 
with the rows ordered in this fashion is scrutinized, it will be noted 
that there is a Aura rearrangement of the columns which wili yield 
a matrix compact along a line--- for square matrix this line is the 
principal diagonal and for the rectangular matrix this line runs from 
corner to corner, The results will be a matrix such as Fígure 9, 

By some criteria such as the type commanders decision that job A 
is most important, the first job order is selected for scheduling, Other 
criteria for finding the starting place are to use the job order row with 
the longest length or the middle row, Having selected the starting row 
by some criteria the complete time scheduling of job order A limits 
certain aspects of job orders B and C (see Figure 9). The time scheduling 
of B and C in turn limit D and E, This continues until a point is reached 
where F and G are the job orders being time scheduled, They can be 
scheduled independently of each other limited only by the job orders 
already scheduled, This is because they use none of the same men, 
machines or facilities. At this point the remainder of the upper half 


can be scheduled without concern for the lower half causing interference, 





Pd p< 
рым 
PS b4 B4 b4 
ма ра ра ра ра 
ра ра ba bed 


<> >< 


P4 b4 bd b< 
ра р В Б 


>< 


u EL en en a Oar та Ф с tr a 


хтазчі рэЗаесівец 6 ә10514 


b< 
b< pd 
Pd bd pq b< bd 
P4 P4 b4 b4 b4 
рс БА рс Bd Pd >d bd 
>< 
t 
і 
i 


ра ра bd be 
>< 


PS 


ра рч рч М М 


> ра bd >d >< 


b< 


p< p< 


SOTIFTTIRH 'seupuow4 ‘ua paz] [ejoeds 


c 


mAM tO x] 


29 


БЈәрлхо qor 








APPENDIX 


The complete prosram is included here for reference purposes. The 
reader 1s also referred to the Flow Charts for individual segments of 
the program, 

The program occupies cells 1200 to 1747 in the main memory. The 
program may be operated either minimum access or sequential numbering 
drum, As written the program requires also the use of a magnetic tape 
having blocks 500 to 620 available in both mode 1 and mode 2, The nonzero 
elements are stored starting in cell 0000 up to 1077, Channel 11 of the 
memory is reserved for temporary storage during the last section of the 
program--- the printout of the correspondence of columns in the undiagonal- 
ized and diagonalized matrices, If part of channel 11 is required for 
the storage of nonzero elements then certain cells and commands must 
be altered for the last segment of the program, (Constant in 1630, 
and command 1664 M. portion). These could be changed to 1200 wiping 
out & portion of the program no longer required, If there are less than 
477 (octal) nonzero elements the program can be simply modified to 
eliminate the use of the tape unit, The partially ordered P's then are 
Stored starting in 0500, The last iteration is then not saved, 

The only constants which must be filled into program are cells; 

1440 --- Number of last cell filled with a nonzero element 


1441 --- Order of matrix plus one, 








If P's are initially completely filled out start computing in 1200, 


If only the Co and Ro numbers are in the м! portion of the P!s start 


computing in 1550, 


lteration is completed, 


completing an iteration, 


AND 


8 1200 
1201 
1202 
1203 
1204 
1205 
1206 


1207 
1210 
1211 
1212 
— 1213 
1214 
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| 
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И 
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A ai e i 


INST 


50 
35 
зЁ 
52 
32 
35 
32 


мі 


С 
144 3 
15506 
2000 
1140 
2000 
2000 


2100 
0000 
1440 
1440 
1612 
2100 
0000 
2006 
ХХХХ 
2007 
2005 
2007 
1222 
1217 
1462 
2006 
1226 
1462 
1226 
3400 
1222 


2007 
2005 
3100 


1460 
145) 
1443 
1444 
1450 
1464 
144% 


2100 
0100 
2100 
1444 
1471 
1615 
2100 
1632 
2100 
1632 
2004 
2100 
1625 
1612 
1217 
2100 
1625 
1612 
1461 
2100 
1615 


2100 
2100 
2100 
2100 


un 


1443 
1440 
2000 
1461 
2100 
2000 
1456 


1477 
0500 
1462 
1226 
1217 
1222 
2006 
2104 
2007 
2105 
123 


1222 
1217 
1217 


1226 
1452 
1215 
1240 
1234 


2006 
2004 
1223 


Ce 
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If SEN 2020 up program will type out all P's after an 


SEN 2040 must be up for program to recycle after 


REMARKS 


CNE TIME CONSTANT PREP 

ln Mt 

N41 in Mi 

Comparison const to 
order 6 S's 

Comparison const for 

tape cmds 


ORDERING 6 S'g 
Initialize partial 2 D? 
Tape block search BL500 
Initialize X 

Initialize put away P! 
Initialize lookup P* 
Initialize put away P* 
Shift P! to 2006 

Extract S! 

Shift P? to 2007 

Extract S? 

Is 52у g1 

Put away P? 

Alter lookup P? 

Alter put away P? 
Tested all P's? 

Put away P? 

Alter put away add P? 
Reduce X to X-1 

Ordered enough S's? 

To next segment of program 
Extract put away add from 1222 


е Я to Pt 
Shift S? to S? 








ADD INST мі мё MS REMARKS 
MAKE HYPOTHESIZED ROW FIRST 
TRANSFORMATION 
1240 32 2100 1444 1241 Initialize Buffer load cmd 
12411 05 3000 3000 ХХХХ Shift partially ordered P's 
1242 15 0000 0101 0000 to tape starting BL500 
1243 35 1241 1467 1241 in Mode 1 
1244 зі 1456 1241 1241 
N 1245 32 1443 1471 1462 Initialize X to N 
Br 21006 32 1462 1471 1250 Prep lookup of P 
E 1247 52 1462 1471 1251 Prep print cmd 
1250 21 9.09 2100 0001 Prints P 
51 %% ХХХ 2100 2000 Shift P to buffer 
1252 32 2000 1450 2107 Extract C to 2007 
1955 32 2000 1551 210h Extract В to 200h 
1254 32 2000 1447 2101 Extract S to 2001 
1255 30 2001 1453 2001 Shift S 4 octal places 
1256 36 2001 1441 1457 Form Y 
(с)—1257 56 2004 1457 2003 Form J-Y 
1260 35 2003 1452 2003 Form K 
1261 36 141 2003 2005 Form K_ 
1262 30 2005 1463 2005 Shift K left 2 octal places 
1265 36 1441 2004 2006 Form J_ i 
1264 50 2006 1463 2006 Shift J left 2 octal places 
1265 32 2100 1444 1312 Initialize put away 
1266 32 2100 1571 1267 Initialize lookup 
1267 35 XXX 2100 2000 Shift P to 2000 
| 1270 32 2000 1445 2102 Extract j to 2002 
' 1271 36 2006 2002 2001 Form (J-4) 
1272 34 2001 2100 1275 Is (J-J) greater than О 
1273 32 2005 1611 2000 If =0, j becomes K 
oh 34 3400 2100 1300 
| 1275 зб 2005 2002 200 Form (K-3) 
| 1276 34 2001 2100 1300- Із (K-j) greater then О 
ZT 32 2006 1611 2000 If -0, j becomes J 
j 1500 32 2000 1614 2102с- Extract 1 to 2002 
| 1301 36 2004 2002 2001 Form_(J-1) 
|! 1302 34 2001 2100 1505 Is (J-1) greater than O 
1303 32 2003 1623 2000 If =0, 1 becomes K 
1304 34 3000 2100 151 
1505 36 2005 2002 2001 | Form (K-1) 
1506 34 2001 2100 1310—— Is (K-i) greater than O 
1507 32 2004 1623 2000 If =0, 1 becomes J 
1310 35 1267 1454 1267< Alter lookup for P 
| 1511 34 3000 2100 11400 To Sum and partial D? S,R. 
| 1312 3 2000 2100 XXXX Put P away 
| 1515 3 1312 1452 1312 Alter put away cmd 
(1314 پاچ‎ 1442 1267 1267 Operated on all P's? 
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ADD INST м1 м M? 
1316 21 1237 1465 0001 
1317 34 1475 1477 E 
1320 34 3000 2100 1421 
1321 05 3000 3000 1500 
1322 30 1457 1463 2006 
1323 56 2007 2006 2006 
1324 55 2006 1464 2006 
1525 50 2006 1466 2005 
1326 56 1441 2005 2005 
1327 50 2007 1466 2003 
1330 56 1441 2003 2003 
1331 32 2100 1444 1356 
1332 32 2100 1471 1333 

[21555 59 XXXX 2100 2000 

| 1334 32 2000 1611 2102 

I 1335 56 2007 2002 2001 

| 1536 зі 2001 2100 1341 

| 1337 32 2006 1616 2000 
1340 34 3000 2100 13444 - 
1341 36 2006 2002 2001 
1342 3% 2001 2100 1344 —- 
1343 32 2007 1616 2000 
1344 32 2000 1625 21020 
1345 36 2005 2002 2001 
1346 34 2001 2100 1351 
1347 32 2005 1626 2000 
1350 3} 3000 2100 155h-—7- 
1351 56 2005 2002 200 
1552 34 2001 2100 1354 
1353 32 2003 1626 2000 
1354 35 1333 1612 133 
1555 34 3000 2100 1400 

| 13556 55 2000 2100 XX 

| 1557 35 1556 1625 1356 
1560 34 1442 1555 1555 
1361 21 1476 2100 0002 
1362 21 1232 1h65 0001 
1565 3 1752 1477 =| 
1364 3, 3000 2100 1421 


REMARKS 


ТОУ ОТ 1476 2100 0002 Print new 3 D* To 
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Tab 
Is previous ZT” greater 
than this X D?, If ote 
go to 1434, 
To subroutine Reload P's 
Return buffer contents saved 


MAKE BYPOTHESIZED COLUMN 
FIRST TRANSFORMATION 
Shift Y 2 octal places 
Form I-Y 

Form L 

Shift_L 2 octal places 
Form L 

Shift I 2 octal places 
Form I 

Initialize put away 
Initialize lookup 

Shift P to 2000 

Extract i to 2002 

Form (I-1) 

Is (I-i) greater than O 
If =0, 1 becomes L 


Form (L-1) 

Is (L-i) greater than O 
If =0, 1 becomes I 
Extract jo to 2002 

Form (1-3) 

Is (1-3) greater than O 
If =0, j becomes L 


Form (L-3) 

Is (L-3) greater then O 

If =0, j becomes m 

Alter lookup for P 

To Sum and partial D? S.R. 

Put P away 

Alter put away cmd 

Operated on all P's? 

Print new x D? 

Carriage return 

Is previous X D? greater 
than this E DS, If greater 
go to 1434 

To subroutine Reload P's 








ADD INST м? Me M9 REMARKS 
1565 3 1457 1452 1557 Make Y=Y-1 
С) 1566 05 3000 3000 1500 Return buffer contents saved 
(8) 1567 21 1237 1465 0001 Tab : ; 
| 1370 34 3000 2100 1571 
1352088 35 1h57 1h52 1257 Ів (Y-1) greater than O 
| 1372 36 1462 1454 1462 Make X-X-1 
| 1373 21 1252 1465 0001 Carriage return. 
| 1374 30 1443 1542 1455 Form N/2 
1375 33 1h62 1455 1216 Is X greater than N/2? 
1576 21 1532 1452 0004 Print MATRIX 1S DIAGONALIZED 
1577 3% 3000 2100 163677 To Correspondence of C orig 
and C new routine 
SUM AND PARTIAL D? SUBROUTINE 
—1400 30 3000 1610 2001 Prepare exit 
1401 32 2001 1615 1415 L i 
1402 32 2000 1614 2101 Extract R 
1403 30 2001 1617 2001 Shift 2 octal places 
1404 32 2000 1616 2102 Extract C 
1405 35 2002 2001 2001 Add C+R 
1406 30 2001 1620 2002 Shift new S LEFT 2 octal 
1407 32 2002 1447 2000 Extract new S into P 
1410 30 2001 1621 2001 Shift 8 RIGHT 2 octal 
1411 36 1441 2001 2001 Form (n+1)-S=D 
1112 26 2001 2001 2001 Form 
1413 35 1477 2001 1477 Add to partial £ D? 
1414 37 1477 3000 1416 Test overflow 
<—1415 3h 3000 2100 XXXX Exit to main program 
1416 35 1476 1452 1476 Tally overflow 
1417 35 2100 1477 1477 Wipe out overflow bit 
1420 34 3000 2100 1415 
RELOAD P's SUBROUTINE 
——1421 30 3000 1446 2001 Prepare exit 
1422 32 2001 1444 1433 Z а 
1423 1 0000 0100 0500 Tape block search BL500 
1424 oh 3000 3000 1500 Save buffer contents 
1h25 32 2100 1hhh 1427 Prepare buffer cmd 
1426 16 0000 0101 0000 Read tape 
1427 ol 3000 3000 XXXX Buffer out to memory 
1h30 55 1427 1467 1127 Add 10 to buffer out addee 
1431 зі, 1456 1427 1426 Have we read gan tape? 
1432 35 2100 2100 1477 Clear partial È 
——>1433 34 3000 2100 XXXX Exit to main program 
——>143h 35 1473 1452 1473 Tally iteration 
1435 35 1476 2100 1,74 Shift partial> D° 
1436 35 1477 2100 1475 Shift partialz D? 
1h37 зі, 3000 2100 1510 








ADD INST mi ve M9 REMARKS 
CONSTANTS 
1440 00 0000 0000 N Number of P's 
1441 00 0000 0000 n+1 Order of matrix +1 
1442 00 N+1 0000 0000 Comparison constant 
1443 00 N 0000 0000 Comparison constant 
1hhh 00 0000 0000 ҮТЇЇ Extractor for M? 
1445 00 0000 0000 7700 Extractor for C 
1446 02 0000 0000 0030 Shift control right 8 octal 
1447 00 0000 0777 0000 Extractor for 5 
1h50 00 0000 0000 7700 Extractor for C 
1451 00 0000 0000 oon Extractor for R 
1h52 10 0000 0000 0001 int control; Constant 1 
1453 02 0000 0000 0014 Shift control right I octal 
1454 10 0001 0000 0000 Print control; Constant 1 
1h55 00 0100 0000 0000 Print control 
1456 00 3000 2000 XXXX Comparison constant 
1457 00 0000 0000 Y Distance between rows (cols) 
being interchanged 
1460 00 0000 0000 0030 Shift control left 8 octal 
1461 00 2006 2100 XXXX Comparison constant 
1462 ОО X 0000 0000 Variable constant 
1h65 00 0000 0000 0006 Shift control left 2 octal 
1464 ОО 0000 0000 0100 
1465 10 0100 0000 0000 Print control 
1466 02 0000 0000 0006 Shift control right 2 octal 
1467 00 0000 0000 0010 
1470 ОО 7700 0000 0000 Extractor for Co 
1h71 00 ТИТ 0000 0000 Extractor for Mi 
1172 00 0077 0000 0000 Extr&ctor for Ro 
1473 оо 0000 0000 0000 Tally of iterations 
1474 00 0000 0000 0000 > of last iteration 
1475 оо 0000 0000 0000 > D? of last iteration 
1476 00 | 0000 0000 0000 > DF of current iteration 
1477 ОО 0000 0000 0000 s. D of current iteration 
1500 00 0000 0000 0000 Buffer save space 
1501 00 0000 0000 0000 " Ё ? 
1502 00 0000 0000 0000 4 у 
1503 00 0000 0000 0000 n 2 i 
150À 00 0000 0000 0000 P E г 
1505 00 0000 0000 0000 n А n 
1506 00 0000 0000 0000 M 2 б 
1507 00 0000 0000 0000 E 7 y 


bs 





REMARKS 


1510 21 1232 1465 0001 Carriage return 


ADD INST M2 М? MS 

1511 21 1473 2100 0003 
1512 21 1232 1465 0001 
1513 17 2020 3000 1525 
1514 1% 0000 0100 0500 
1515 38 2100 1444 1516 
1516 05 3000 3000 XXXX 
DONT 05 0000 0102 0000 
1520 35 1516 1467 1516 
1521 34 1456 1516 1516 
1522 17 2040 3000 1207 
1523 22 0000 0000 0000 
152% 3% 3000 2100 1200 
1525 35 1440 1452 2001 
1526 32 2001 144) 1527 
1527 21 0000 2100 XXXX 
1530 21 1232 1465 0001 
15531 5h 3000 2100 1514 
1532 00 3476 4175 6360 
1533 00 5664 6044 6447 
1534 Оо 6141 7377 6641 
1535 00 7160 5545 4731 
1536 35 1176 1452 1476 
1537 35 2100 1477 1477 
1510 5% 3000 2100 1575 
1541 00 0000 0000 0000 
1542 02 0000 0000 0001 
1543 00 0000 0000 0000 
1544 00 1633 0000 0000 
1545 00 0000 0000 0000 
1546 00 0000 0000 0000 
1547 00 0000 0000 0000 


Print iteration and p? 

Carriage return 

SEN 2090 UP, print P's 

Tape block search BL50O 

Initialize buffer cmd 

Buffer load 

Tape write 

Add 10 to buffer cmd 

Have we written enough tape? 

SEN 2040 UP, recycle 

Halt 

After halt, will recycle on 
compute command 

Form N«1 

Extract (N41) into print cmd 

PrintP's 

Carriage return 


Flexo code 
$” п 


Tally an overflow 
Wipe out overflow bit 


Buffer save space-minaccess 
Shift control right 1 binary 
Buffer save space-minaccess 
Constent 

Buffer save space-minaccess 


Buffer save space-minaccess 


Only the M portion, 1.e,, Co and Ro, need be initially put in 


the computer, If only the мі portion is filled, starting computing 


in cell 1550 will complete the other portions, namely, current R and 


C and S for all elements P, 


compute ar initial EXE If SEN 2010 up thisZ D 


This portion of the program will also 


^ ‘111 bel ESS 


and the initially filled out P's saved in tape unit. 








ADD INST M1 M? M3 REMARKS 
COMPLETION OF INFORMATION 
FOR P's: COMPUTE INITIAL? D*? 
1550 55 2100 2100 1477 Initialize ial 5 р? 
1893 30 1440 1560 1443 Form N in M 
1552 35 1443 1454 1440 Form N41 in M 
1553 52 2100 1571 1555 Initialize lookup cmd 
1554 32 2100 1444 1575 Initialize put away cmd 
1555 55 ХХХХ 2100 2000 Shift P to buffer 
1556 52 2000 1470 2102 Extract Co 
1557 52 2000 1472 2101 Extract Ro 
1560 50 2001 1546 2001 Shift Ro 8 octal places 
1561 50 2002 1446 2002 Shift Co 8 octal places 
1562 35 2002 2001 2003 \ Extract С апа R 
1563 32 2003 1444 2000 ) into P 
1564 30 2001 1463 2001 Shift R 2 octal places 
1565 35 2002 2001 2001 Form C+R 
1566 30 2001 1463 2002 Shift S 2 octal places 
1567 32 2002 1447 2000 Extract S into P 
1570 30 2001 1466 2001 Shift S right 2 octel 
1571 56 2001 1441 2001 Form D=S-(n+1) 
1572 26 2001 2001 2001 Form 
73 55 1477 2001 1477 Add D? to partial > D? 
1574 37 1477 2000 1536 Test overflow 
1575 35 2000 2100 XXXX Put away P- 
1576 35 1555 1454 1555 Alter lookup 
ET 55 1575 1452 1575 Alter put away 
1600 34 1442 1555 1555 Have we operated on all P? 
1601 21 1476 2100 0002 Print = 
1602 32 1440 1450 2100 Make comparison constant 
1603 35 2000 1464 2000 for buffer-tape cmd 
1604 32 2000 1444 1456 
1605 17 2010 5000 1434 SEN 2010 UP, shift D? and 
store P's. DO NOT HAVE 
SEN 2040 UP 
1606 22 0000 0000 0000 
CONSTANTS OPTIMUM SPACED 
1607 00 0000 0000 0000 
1610 02 0000 0000 0030 Shift control right 8 octal 
1611 00 0000 0000 7700 Extractor for C 
1612 00 0001 0000 0000 Print control; const 1 
1613 00 0000 0000 0002 
1614 ОО 0000 0000 0077 Extractor for R 
1615 00 0000 0000 7777 Extractor for M? 
1616 00 0000 0000 7700 Extractor for C 
1617 00 0000 0000 0006 Shift control left 2 octal 
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FLOW CHART - CORRESPONDENCE OF ORIGINAL AND NEW COLUMNS 


[Initialize LM ii 


La] 


Shift P to buffer, extract S . 
(Form S-(n+l) =D — 


m uoo valeo ep cmm дй 





-- -—— 5 4 -— 5. c Pme cue ME, шама ———. S m NS 


un Ga San u en u er = а 8 ә з 


Put P. P away 7 in temporary storage | 
Alter put away to. (а) 





.. po . e... 5 
evt * 
per? ) 
Cm. Pun JB 


——Yege-«/ TM2O » ———3 No 
E p 


Have we tested all P 
i.e0., N15» lookup address 





$—— — ee te to em à 


(Form extractor of C c orig. | 


з шы Зь ль O rn BE ne dern Ф ў уў ta Кї, A a 


Initialize Put away Dy to (a+n) 
Form ee. — 


A PETALA = 





TE gel ee ge, 


Initialize lookup .P^, put avay P3 | 


pou. we 
| Shift P! to buffer, extract C 
Shift P? to buffer, extract C? 
No &-———— С^» C. Yes E s 
Put pP? in | Put E! in pa “storage | 
P2 storage Shift P? to P', 252 to St 
ша . .. .. ipic M e BEEN AN 


‚Alter lookup of P^ 1 
¡Alter put away of р я 


Af. SOMIT. еа дага 





A а-а a нда 





Have we tested all P's 
(a+n\>lookup 











Yes « 


Put PT in là3st storage position  . 
Alter last storage cel! to (a n)-l 


No ds cina Ф‹ (ы ' 4” э, тые бай 14% Wir ee ست کے سیک ب‎ &—‹“%Р ъ 


—-—1s last storage cell 2 (a+l) 


2 ar EEE gr Sg uy eray naan A, > a m 
| ds =: y t m ts OM gt er аачы MED афа. 


a паь Чаран A A A с.-х ыё en AP, Pt alis ^ a aE ——- — _ ————  — a —_—-— —— 00 


gma hee _ 
* 


©4, 4 wee onu rro m -— 


Form extractor of C new | 
Initialize lookup | 
*Print heading "Corig-Cnew" | 


A — © Were, af ante a 


De DA Pr pm D yde W کبک‎ 





Shift P to buffer cell 2000 
Extract C, Shift left € octal 
Put in cell 2001 
Print 2 words-first 4 characters only 
Alter lookup | 


Maa? ru, EP tb oue Re 


- 4 ue P amos che acum He a سے چس سے یویم فی‎ че ць атавы Ф Ф цы nego e ha PA O ULA вы TEA o Ge cm Ios ga ERR ARP UR A 





“я 


——— Is (am)>Address of P? 


| * 
[Add 252525252525 to itself |j 
hoor 
юе Olo — -o EE 
Y 





f њо За (р чо отеу о 


Alter cmds marked * 
, Print Headings | 
| "Cnew-Corig' | 
| Print conirol to і 
И print 2007 and 2000 
| 

| d 





Shifted C is put in 
2007 


A c E‏ و س ل 


е 


Figure 10 


З= 


¡Alter cmds marked * ` 
back to original 
ı settings __ | 


rr A o s e Rng + © a Tg 


Aue o e mm de A arte i 


| 


| سے ا | 


¡End of Program __ і 


سید 





1660 


T 


4 


з558888535395358 


мі 


0000 
0000 
0000 
0000 


0000 
0000 
2525 
1100 
1736 
0000 
3436 
6360 
6232 


11441 
1447 
2100 
1630 


2000 
2003 
2001 
2000 


1646 
1642 
1442 
1470 
1630 
2000 
2001 
2001 
2001 
1630 


Me 


0000 
0777 
0000 
0000 
0000 
0000 
0000 
2525 
0000 
0000 


6232 
7327 
6466 


1622 
2100 
1471 
1444 
2100 
2004 
2007 
2100 
2100 


1452 
1454 
1642 
2100 
1941 
1615 
1444 
1452 
1460 
1444 


39 


REMARKS 


Shift control left 2 octal 
Extractor for S 
Shift control left 4 octal 
Extractor for R 
Shift control right 2 octal 


Extractor for R 
Test constant 
Constant for temp storage 
Constant 
Extractor for 8 
Flexo code 
11 ?f 


11 " 


Form (n+1) in M° 

Shift extractor of C to buffer 

Initialize lookup 

Initialize put away 

Shift P to buffer 

Extract S into 2003 

Form S-(n+1)=D 

Is D greater than O 

If =0, P is a diagonal element 
put in temp storage 

Alter put away 

Alter lookup 

Have we tested all P's? 

Shift extractor of Co to 2002 

Form 1100+(n+1) 

Form 1100+(n-1) 

Initialize put away P! 

Form 1100+n 

Shift left 8 octal 

Make & comparison constant 





1661 


—> 1662 


1663 
155% 
1665 


— 1446 


1667 
1570 
167] 


1672 





== 


IO ےھ‎ аа = 


1673 
1574 
175 
1676 
1677 
1700 
1701 
1702 
1703 
1704 
1705 
1706 
1707 
1710 
1711 
1712 
1713 


r>171% 


1715 
1716 
1717 
1720 


1721 
1722 


1723 
1724 
1725 
1726 
1727 
1750 
1751 
1732 
1733 
1734 
1755 
3736 
ET 
1740 


1630 
2000 
1630 
1100 
2006 


2007 
2005 
2007 
1666 
1671 
2003 
2006 
1675 
2003 
1675 
3000 
1671 
2006 
2007 
2005 
3000 
2001 
1633 
1232 
1450 
1630 


2000 
2001 
2000 
1714 
1232 
2003 
1627 
1627 
1631 
1220 
121] 
5000 
2100 
1544 
1252 


0000 
3436 
4551 
6477 


е 


145% 
1471 


M? 


2000 
1666 
1671 
2006 
210h 
2007 
2105 


17028 


1666 : 


1 
$ 
\ 


! 


162] 


1666 


1675 | 
2003 


1662 


17077 


` 


1703 


2006 
2004 
1672 


2005«— 


0003 
0001 
2002 
1714 
2000 
2101 
2001 
0002 
1714 
0001 
1714 
1627 
1731 
1710 
1717 
1716 
165% 
1627 
1710 
1717 
1716 
0000 
6466 
6232 
(EN 


FEMARKS 


Make & comparison constant 
Initialize lookup 
Initialize put away p? 
Shift Pl! to buffer 
ExtraOt C! into 2004 
Shift P? to buffer 
Extract C? into 2005 

Те Os cn 

Put away P? 

Alter lookup P? 

Alter put away P? 

Tested all P's? 

Put away P! 

Alter put away Р! 

Alter comparifon constant 
A11 P's ordered? 


Prepare put away P! 
Put Pi away 

Shift P? to P! 
Shift S? to 51 


Prepare a comparison constant 
Print headings 

Carriage return 

Shift extractor for C to 2002 
Initialize lookup 

Shift P to buffer 

Extract C 

Shift C left 8 octal 

Print Co =- C 

Alter lookup 

Carriage return 

Have all C's been printed? 
Make tally constant 

Is this 2nd time thru routine? 
Change print control 

Change print control 

Change C put away address 


Wipe out overflow bit 

Change print control 

Change print control 

Change C put away address 

Final halt 

Flexo code 
п 
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